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Abstract 

The braiding of the solar coronal magnetic field via photospheric motions - with subsequent 
relaxation and magnetic reconnection - is one of the most widely debated ideas of solar physics. We 
readdress the theory in the light of developments in three-dimensional magnetic reconnection theory. 
It is known that the integrated parallel electric field along field lines is the key quantity determining 
the rate of reconnection, in contrast with the two-dimensional case where the electric field itself is 
the important quantity. We demonstrate that this difference becomes crucial for sufficiently complex 
magnetic field structures. 

A numerical method is used to relax a braided magnetic field to an ideal force-free equilibrium; that 
equilibrium is found to be smooth, with only large-scale current structures. However, the equilibrium 
is shown to have a highly filamentary integrated parallel current structure with extremely short 
length-scales. An analytical model is developed to show that, in a coronal situation, the length 
scales associated with the integrated parallel current structures will rapidly decrease with increasing 
complexity, or degree of braiding, of the magnetic field. Analysis shows the decrease in these length 
scales will, for any finite resistivity, eventually become inconsistent with the stability of a force-free 
field. Thus the inevitable consequence of the magnetic braiding process is shown to be a loss of 
equilibrium of the coronal field, probably via magnetic reconnection events. 

Keywords: MHD; magnetic fields; Sun: corona. 
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1 Introduction 



The heating of the solar corona as a result of the braiding of its constituent magnetic loops is one of the 
fundamental theories proposed to explain coronal heating. Despite its long history relative to the subject 
area - originating with the early notion of 'topological dissipation' (Parker, 1972) - the braiding concept 
is highly controversial, with no general consensus reached as to the exact mechanism or even its general 
viability. 

Coronal loops are typically modelled as consisting of ideal plasma threaded by a force-free magnetic 
field. With the ends of each loop anchored in the turbulent photosphere and Alfven travel times along 
a loop much faster than the time-scales associated with convection, the field evolution should be via a 
sequence of force-free states. Parker (1972) argues that relaxation to a smooth equilibrium following 
an arbitrary perturbation is impossible (excepting artificial cases with certain symmetries) and instead 
tangential discontinuities, corresponding to current sheets, must develop in the field. The claim was first 
refuted in an analytical treatment by van Ballegooijen (1985) who demonstrated the existence of smooth 
equilibria in non-symmetric 3D configurations and under arbitrary continuous boundary perturbations. 
He found singular current sheets appearing only with the application of discontinuous velocity fields at 
the boundary or in the case of a discontinuous magnetic field, corresponding to isolated magnetic flux 
patches in the photosphere (an idea followed-up in Priest et al. 2002) . 

Several other authors have also demonstrated classes of smooth equilibria (Zweibel & Li 1987, Long- 
cope & Strauss 1994, Craig & Sneyd 2005) and numerous authors have tackled the problem via numerical 
simulations (e.g. Mikic et al. 1989; Longcope & Sudan 1994; Galsgaard & Nordlund 1996) finding equi- 
libria that are non-singular - although they do often contain current structures on a smaller scale than 
that of the footpoint displacements. However, it can be argued that numerical simulations may have 
an inherent resolution problem in allowing current sheet formation to be recognized. In their analytical 
treatment Craig & Sneyd (2005) demonstrated the existence of smooth equilibria for any footpoint distur- 
bance and for arbitrary compressibility Significantly, their argument was further advanced through the 
use of a Lagrangian relaxation scheme that enabled the stability of smooth solutions to be demonstrated, 
even for significant footpoint displacements. 

While no firm conclusion has been reached, either analytically or numerically, on the capacity of 
magnetic braiding to produce current sheets in the corona, it is clear that the present thinking regard- 
ing magnetic braiding is fundamentally based on ideas developed from classical two-dimensional (2D) 
reconnection theory (Sweet 1958; Parker 1957; Petschek 1964) namely that, in order to allow for rapid 
reconnection, extremely small-scales in the electric current (compared with the large-scale magnetic field 
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Figure 1: (left) Illustration of the pigtail braid, on which the braided magnetic field modelled here is 
based. The pigtail braid is the braid equivalent of the Borromean rings (right) 

structure) are required. In three dimensions (3D) a necessary condition for the occurrence of reconnection 
is a non-zero integrated electric field parallel to the magnetic field (Hesse & Schindlcr 1988, Schindler et 
al. 1988). In particular, reconnection may take place in the absence of any topological features such as 
magnetic null points. As such it has quite different characteristics, both in terms of onset and evolution, 
to the two-dimensional case. For some examples of 3D reconnection in the absence of a null point, the 
situation relevant to magnetic braiding see, for example, Hornig & Priest (2003) and Pontin et al. (2005). 
Hence the results of magnetic braiding should be reassessed with the proper 3D reconnection criteria in 
mind. In particular the development of parallel electric fields as a result of braiding processes should be 
examined, together with their consequences. We begin to address these issues in this paper, within the 
framework of resistive MHD. 

The outline of the paper is as follows. We build an analytical model for a braided magnetic field 



between two parallel plates (Section 2.1) and describe a numerical relaxation scheme in which we relax 



this field to a force-free state, as expected for the solar corona (Section 2.2 1. We then analyse the result of 
the relaxation process in terms of both the current structure and integrated parallel currents (Section |3| 
and discuss how these structures depend on the degree of braiding of the field (Section [3} . Results are 
discussed in Section [5] and conclusions drawn in Section [6j 
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2 Methods 



In this paper we are interested in how braided magnetic fields relax to a force- free equilibrium and specif- 
ically in the nature of such equilibria. We begin by describing a general method for constructing braided 
fields and then concentrate on one particular field structure, a configuration based on the Borromean 
rings. This situation has been chosen since it represents a truly braided field with a number of advan- 



tageous solar-like properties (as described in Section 2.1|, lends itself to analytical modelling and, in 
addition, can be modified in a way to allow for easy comparison with more and less complex fields. We 
do not concentrate on the details of how this particular field would be built up in the corona, rather 
taking it as a fairly representative example of a complex field. The construction method for this and 
similar fields is described in the following section. 



2.1 Braid construction 

The braid considered here is modelled on a standard 'pigtail' braid. The pigtail braid consists of three 
strands, braided as shown in Figure [I] (left). By way of motivation for this choice we suppose that the 
braiding of field lines anchored in the photosphere occurs in a random manner, that is, left-handed and 
right-handed twists occur with the same probability. We therefore wish to take a braid which has no net 
twist. More precisely, if we close the braid to form a link, by identifying the upper and lower planes in 
Figure [T] then the pairwise linking of strands should vanish. The simplest link with this property is the 
Borromean rings (see right hand image of Figure [lj which have the property that any two of the three 
rings are unlinked whilst the whole link itself is non-trivial, i.e. cannot be taken apart. Our pigtail braid 
is the braid equivalent of the Borromean rings. The vanishing net twist of our configuration is also the 
most conservative case. It is easy to drive currents in a model by highly twisting the field, a situation that 
is not observed in the corona, whilst our choice represents an essentially twist-free situation. A further 
advantage of the braid is that it can be considered as the concatenation of three identical parts, each 
consisting of a pair twists (intersection of strands), one positive and one negative. As such it is easy to 
study the effect of increasing braid complexity through a concatenation of increasing numbers of these 
elementary parts. 

The magnetic fields representing our braids are constructed as follows. Each field consists initially of 
a uniform background field (b$i), superimposed on which are a number of isolated magnetic flux rings, 
corresponding to the crossings of strands in our braid. The basic form of these flux rings is given, in 
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cylindrical polar coordinates, by 

T ( \ 

B c = 26 fc-exp I -— -— I e , (1) 
a \ cl L j 

with the parameters included allowing for variation in the field strength (k, b ), as well as the radius (a) 
and vertical extent (I) of the region of toroidal field. The effect of adding the field given by equation (|T|) to 
the background field &oz is to create a region of non-uniform twist within the otherwise straight vertical 
field, with the maximum degree of twist controlled by the parameter k. An expression corresponding to 
equation ([l]) may be found in Cartesian coordinates with the field centre transposed from r = 0, z = 
to x — x c , y — y c , z — z c . We label this field B c : 

Tj o b k { ~(x-x c ) 2 -(y-y c ) 2 (z-z c ) 2 \ 

B c = 2 exp — (- (y - y c ) x+ [x - x c ) y , (2) 

a \ a 2 I J 

and let B c . denote the field B c taken with the parameters = (x c ,i, y c ^, z c> i, ki,cii,li). Braided fields 
of various complexity may then constructed as 

n 

B = 6 z + J]B Ci , (3) 

i=i 

with appropriately chosen parameter sets Ci (i = 1 . . .n). In this paper we set bo = 1 throughout. 

The braided fields constructed in this manner are taken as an initial condition in numerical relaxation 
experiments. An alternative would have been to begin with a uniform magnetic field between lower 
and upper boundaries and use a numerical technique to simulate photospheric motions and obtain a 
braided field and to then proceed to relaxation techniques. However, here our primary interest lies in 
the relaxation process and its consequences and so we take the analytical models as a proxy for a field 
constructed through simulated photospheric motions The advantages of this approach lie in a significant 
saving in computational time, the ability to easily adopt fields with the aforementioned braid properties 
and a model that lends itself to analytical analysis. 

To construct the pigtail braid on which we base the majority of our analysis (Sections [2] and [3]) we 
take the parameter sets c x = (1,0,-20,1,^2,2), c 2 = (-1, 0, -12, -1, V2, 2), c 3 = (1,0,-4,1,^2,2), 
c 4 = (-1,0,4,-1,^2,2), c 5 = (1,0,12,1,^,2), c 6 = (-1, 0, 20, -1, y/2, 2) , and let 

6 

B = fo z + ^B Ci . (4) 

i=i 
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(a) Some particular 
field lines within E 3 . 




(b) Illustrating the same field lines as in (a) with a 
different aspect ratio and view-point to show more 
clearly the braided nature of the field. 




(c) Isosurfaces of current, j = j'max/4. 

Figure 2: Initial state for the braided magnetic field E 3 . 
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As previously mentioned the field is the concatenation of three identical parts. The 'elementary' unit, 
consisting of one positive and one negative twist, we label E for convenience and so refer to the pigtail 
braid as E 3 . In a similar manner we may construct the fields E 2 , E 4 etc. The construction of E 3 is 
such that the field lines threading the regions of maximum twist form a Borromean braid as shown in 
Figure |2~T) However only a small proportion of the total flux is braided in this manner; field lines typically 
pass through fewer of the regions of toroidal field, turning by a smaller angle as they do so. The regions 



of twisted field result in six isolated regions of current, as shown in Figure 2.1 ^c). 

The field is defined in infinite space, while to realistically model coronal loops we should take lower and 
upper boundaries, representing the photosphere, to which the field is perpendicular. For the analytical 
work that follows, the lower boundary can be conveniently considered by taking the limit z — > — oo and, 
similarly, the upper boundary as z — > oo. Practically however, due to the exponential decay of the toroidal 
field components, a short distance above and below the regions of toroidal field centered at z = ±20, the 
field can be approximated as vertical. Accordingly, for the numerical investigations we define the domain 
as x £ [— 4, 4],y <G [—4, 4], z E [—24,24]. Taking this domain it is worthwhile to note that the difference 
in magnetic energy between E 3 , as defined by Equation Q, and the potential field satisfying the same 
boundary conditions (i.e. the uniform field b^z) is just 3.08%; the field is close to potential in this respect. 
This is in accordance with observations of coronal fields (Solanki et al. 2006) . 

We now proceed to describe the numerical method we will use to relax E 3 to a force-free state whilst 
keeping the field on the boundaries fixed. 



2.2 Numerical methods 

Whilst, as previously mentioned, the field in the solar corona is expected to be largely force-free, i.e. in 
an equilibrium state satisfying (V x B) x B = 0, the initial state for E 3 , as given by Equation Q is 
not in such an equilibrium. We therefore seek to relax that field towards a force-free equilibrium while 
preserving the topology of the field, i.e. the relaxation must correspond to an ideal plasma evolution within 
the volume and have vanishing velocity on the boundaries. For this we use the 3D Lagrangian magneto- 
frictional ideal relaxation scheme described in Craig & Sneyd (1986). The advantage of a Lagrangian 
scheme is that since the numerical grid and, due to the ideal evolution, the magnetic field lines, moves 
with the plasma elements, grid points accumulate in regions where length scales become small. This 
enables a better resolution of these small-scale features than would be achieved with a similar Eulerian 
code and the same number of grid points. 

Details of the code can be found in Craig & Sneyd (1986, 1990) and are, therefore, only briefly 
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summarised here for convenience. Since the relaxation is ideal, both the line element <5x that joins two 
fluid particles and the vector B/p must follow the same time evolution, that is 



and similarly for 5x; the code makes use of this by invoking the same Lagrangian description for both. 
In order to realize the equilibria, a fictitious equation of motion that guarantees a monotonic decrease in 
the magnetic energy of the system is adopted. The momentum equation takes the form 



i.e. a frictional term proportional to the velocity of the fluid, v, is taken (y ~ 1) and inertial effects 
ignored. The choice is also advantageous in that it is a parabolic system, so the possibility of shock 
waves, and is justified since our primary interest is in the final stable equilibrium rather than the path 
to it. Additionally, this choice of equation has been shown not to alter the stability or state of the final 
equilibrium. Magnetic flux and V • B are both conserved by the relaxation. The numerical technique 
itself is an implicit (ADI) unconditionally stable scheme. 

3 Force-Free State 

We now take the initial state for E 3 (as given by Equation [4]) and use the numerical method described 
in Section |2.2| to relax it towards a force-free equilibrium, with the field on the upper and lower bound- 
aries [z = ±24) held fixed by line-tying at those boundaries. The experiment has been run at various 
resolutions with the results presented here being from a 161 3 computational grid and a relaxed state 
with (j x B) max k2x 10~ 2 (a figure which may be compared with the equivalent initial value of 1.38). 
Numerical accuracy issues prevent further relaxation but we view this state as a good approximation to 
force-free; for a discussion see Section [6] 

The final equilibrium is found to be smooth in accordance with the predictions of van Ballegooijen 
(1985) and others, for example, Craig & Sneyd (2005). Figure [3]ja) shows an isosurface of current, 
in the force-free state. Two spiraling current structures are seen to extend vertically throughout the 
domain. No current sheets (singular or non-singular) are formed in the relaxation, indeed whilst initially 
the maximum absolute value of the electric current is |j| max ~ 2.83, in the final state this value is 
lilmost ~ 1-47. Considering the field structure itself, some field lines illustrative of E 3 in both the initial 




"v = jxB, 
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(a) (b) 

Figure 3: (a) Isosurfaces of current density, (|j| = 0.45) in the relaxed equilibrium, (b) Particular field 
lines in E 3 in the initial state (black) and the equilibrium state (red). 

and final states are shown in Figure ^b) . The same three field lines are plotted in both cases - the 
boundary conditions allow us to identify these exactly. As illustrated in the figure, the twist in the initial 
configuration has become distributed evenly along field lines although they are, of course, still braided 
since the relaxation was ideal. The magnetic energy in the relaxed, equilibrium state is just 0.96% in 
excess of potential, a figure which may be compared with the 3.08% excess of the initial state (as discussed 



in Section 2.1l 



That a smooth force-free equilibrium has been found for E 3 adds weight to the body of evidence 
against the Parker hypothesis (Parker, 1972) and lends support to the arguments of van Balleooijen (1985) 
and Craig and Syned (2005) that smooth coronal equilibria may be found after arbitrary photospheric 
evolutions. However, as discussed in Section [l] in three dimensions it is the integrated electric field 
component parallel to the magnetic field that is the crucial quantity for reconnection. Since within the 
framework of resistive MHD, as considered here, Ohm's law is given by E + v x B = 77 j , the parallel 
electric fields are directly related to parallel electric currents by the relation J Eudl = J rjjudl = rj j j\\dl, 
where the last equality holds in the case of uniform resistivity 77 (in practice, 77 may be spatially varying, 
dependent, for example, on the current density itself). These considerations motivate us to consider the 
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Figure 4: Contours of \J j\\dl\ for E 3 in the initial (left), and equilibrium (right), states. The integral 
is taken along field-lines starting at the given location on the lower boundary through to the upper 
boundary. 

nature of J j\\dl. Results are detailed in this section and their implications for the corona discussed in 
Section [3 

We calculate the integrated parallel currents, f j\\dl, along field lines in E 3 , both for the equilibrium 
state and, for purposes of comparison, the initial state. Specifically, for a large number of points on the 
lower boundary of the domain (z — —24) we integrate the parallel current along the field line starting 
from that location until it reaches the upper boundary (z = +24). Figure [1] shows a contour plot of the 
result for a subsection of the domain. In the figure the contour value indicates the total absolute value 
of integrated parallel current, \J j\\dl\ , for the field line starting at that given location. Two important 
points about the structure of the integrated parallel currents are illustrated: 

1. The width of the layers of enhanced \ f j\\dl\ is very small compared with both the size of the domain 
and the typical scales of j and B. This statement is made more precise in the following text. 

2. The structure, or distribution, of JjmcU is approximately preserved in the relaxation process, 
although its peak value increases slightly (from 27.04 to 33.86). 

The second of these points, the near preservation of \J j\\dl\ along field lines during the relaxation 
process, has the consequence that the key properties of the equilibrium field may be deduced from the 
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initial, non-equilibrium state. This allows us to examine the distribution of \J j\\dl\ in various braided field 
configurations using only the initial state. This state is imposed as a (physically motivated) analytical 
dependence and so we may quickly and accurately determine the distribution on the lower boundary of 
integrated parallel currents along field lines without resorting first to computationally expensive relaxation 
experiments. 

We focus now on the second of our key points, the width of the integrated parallel currents. Figure [5] 
gives contours of that quantity in the initial state (comparable to the equilibrium state as discussed 
above) for both the full domain and for a small subsection of it. The integrated parallel currents are 
typically concentrated in thin layer-like structures, the exception being the two broader peaks evident 
in the left hand image of Figure [5] The right hand image of this figure shows one such typical layer in 
some detail. This layer has a half-width at half-maximum of 0.0167. The width be compared with the 
typical scales of the electric current, j, in the equilibrium - in the horizontal direction the half- width at 
half-maximum of the current layer is ~ 0.74; the typical scale of the integrated parallel currents is 3% of 
this value. We address the possible consequences of the existence of these small scales in the integrated 
parallel current in Section [5] 

Recall that in resistive MHD the parallel electric fields and parallel electric currents are related via 
J E\\dl = Jr)j\\dl, and that, in the case of uniform resistivity, J j\\dl (as considered here) will give us the 
structure of J Ei\dl. 

In practice, the resistivity, 77, may not be spatially uniform but dependent, for example, on the current 
density, j. In the force-free equilibrium found we have |j| ~ \B\ and, in addition, \B\ is dominated by the 
constant B z component. Hence the current here is essentially constant along field lines. In the perhaps 
more realistic case of a non-uniform resistivity, a dependence such as rj ~ \j\ n (n > 0), for example, would 
therefore lead to integrals like J j n+1 dl. These would show even more enhanced structures and so taking 
77 as constant we are considering the most conservative case. 

We now proceed to consider how changing the the complexity of the braided field affects results. 

4 Effect of Braid Complexity 

The initial state in E 3 , as given by Equation [4] consists of six isolated magnetic flux rings (given in that 
expression by B c . , i = 1 ... 6) superimposed onto a uniform background field (&oz), the effect being to 
create six regions of isolated twist within that otherwise straight field. As discussed in Section [2TT] fields 
similar to E 3 may also be defined, using Equation ^ and the appropriate choice of parameters (guided 
by Equation [4]) to maintain the same essential pattern in the braid, including the alternating sign of 
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Figure 5: Contours indicating the value of | / j||d/| along field-lines in the initial state of E 3 for (left) the 
full domain and (right) a small subsection only. 

twist with height, but with varying numbers of twisted field regions. More precisely, various fields E l 
(i = 1,2,3,...) may be constructed by reference to the 'elementary' unit, E, with magnetic field given by 

2 

B = 6 z + 5Z B6 " 

i=i 

where Ci = (l, 0, — 4, 1, v2, 2) and &2 = (— 1, 0, 4, — 1, v2, 2), and choosing a suitable composition. We 
name these fields according to the number of concatenations taken, so E 2 consists of four twisted regions, 
E 4 of eight, etc. 

We now analyse the pattern of integrated parallel currents along field lines for a number of fields 
E % in a similar manner to that carried out for E 3 . The aim is to deduce how the width of the \J j\\dl\ 
structures varies with braid complexity. Our field construction method allows the complexity of braided 
field E n to be denoted by n; we take an increase in n to be generally representative of an increase in the 
degree of braiding of the coronal field due to additional photospheric motions. 

Figure [6] shows contours of \J j\\dl\ for portions of the domain for both E and E 5 . The scale of the 
domain shown in the pictures has been chosen to allow the length-scales of | J j\\dl\ to be clearly seen. For 
E the maximum value of \ J j\\dl\ is lower (\J j\\dl\ = 11.38) Additionally, the characteristic widths 
of the distribution are much broader. For E 5 the peaks in |J i||d/| are much narrower (as evidenced by 
the scale of the domain shown) and have a higher value (\J j'||d/| = 33.62 in the section of the domain 
shown which has been arbitrarily chosen and, as such, is unlikely to encompass the global maximum in 
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Figure 6: Contours of |/j|idi| along field lines for appropriate portions of the domain in (a) E, and (b) 




n 



Figure 7: Scatter plot showing the variation in the half-width at half-maximum (d) of the peak in \f j\\dl\ 
along the line y = x (0 < x < 2) with n, where n denotes the field E n . The width of the integrated 
parallel current structures is seen to decrease exponentially with braid complexity - the line is given by 
d = k1CT™ /2 , where k = 0.77. 
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From these simple considerations we expect that as the complexity of a braided field increases then 
the maximum in J" j\\dl\ increases while the width of the J j\\dl\ structures, d, decreases. 

In order to determine the functional dependency d(n) and the maximum value of \f j\\dl\ with n 
(number of elementary braids, E) we would have to make two-dimensional contour plots such as those 
in Figure [6] for many higher numbers, n. Due to the increasing resolution required for these plots such a 
method is numerically very costly. Instead we choose to consider integrated parallel currents along field 
lines originating along the line y — x, < i < 2 on the lower boundary. Taking just a ID line cut has 
the disadvantages that the line may not be (indeed is unlikely to be) perpendicular to \J j\\dl\ structures 
and may not cross regions of maximal \J j\\dl\. However the method might still be sufficient to allow a 
rough functional dependence of d(n) to be determined. For a number E n we consider the half-width at 
half-maximum, d, of the maximum peak in the integrated parallel current value along this line. We note 
that the value obtained may not be representative of the thinnest structures within the full domain but 
the method should suffice in giving us an indication of the nature of the variation of d(n) . 

The result is shown as a scatter plot in Figure [7J we deduce that the width of the integrated parallel 
current structures decreases exponentially with increasing braid complexity - the line superimposed is a 
fit to the data with equation d = K1CP™/ 2 , where k = 0.77. Although this equation should not be taken 
to given any exact dependence of d with n, it gives a rough guide as to what we may expect. In addition, 
we note that the distance between adjacent peaks in \ J j\\dl\ also becomes smaller with increasing n. 

The above described dependence of d with n might be expected from the following theoretical con- 
sideration. If we assume that the maximum value of \J j\\dl\ in E coincides with the location where the 
field lines are stretched the most then we would expect d ~ do/X where A is the eigenvalue of the field 
line mapping and do the length-scale of the current distribution. Subsequently an iteration of this map, 
as appropriate for a concatenation of E to form various braided fields E n , would result in the width d 
of \ J j\\dl\ layers in E n of d ~ do/\ n . This would occur provided there are sufficiently many regions of 
maximal field line stretching that it becomes likely that a region with this property is mapped onto a 
similar region in the appended braid. By the same argument we would expect the maximum value of 
| J" j\\dl\ to increase linearly with n. 

The implications of the presence of narrow \J j\\dl\ layers within a magnetic field structure are dis- 
cussed in the following section. 
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5 Implications for Coronal Evolution 

Our results show clearly that small-scale electric current structures are not a necessary requirement 
for small-scale integrated parallel current structures. The consequences of such a situation are less 
clear. Whilst is is well known that in two dimensions thin current sheets can lead to rapid magnetic 
reconnection, an equivalent condition for the three-dimensional case has not been deduced as yet. In this 
section we show that small scales in the integrated parallel currents are a likely candidate; they are, at 
least, inconsistent with the presence of a stable force-free non-equilibrium. 

We used a numerical code to relax E 3 towards a force-free equilibrium state, j x B = 0. In such an 
equilibrium the ideal Ohm's law, E + vxB = 0, may be satisfied with v = and 4> — 0, where E = — V</> 
due to the stationarity of the state, a stable force-free ideal equilibrium. In a real resistive plasma, such 
as the solar corona, a very small but finite resistivity, r), is present and Ohm's law may be given by 

E + vxB = 77j. (5) 

Then, since the equilibrium is time independent, E = — V0 and the electric potential <f> may be deduced 
by integrating along field lines. For a spatially uniform 77, <f) is given by 

4> = ~ J V3\\ dl = -v J dl. 
The plasma velocity v then follows as 

(— V0 - r/j) x B (-V0) x B 



v 



B 2 B 2 



where the second equality holds since the equilibrium field is force-free. From this we estimate the plasma 
velocity as 

since the narrow width d of the parallel current structures will induce the largest gradients in the electric 
potential. Similarly, making an estimate of the size of the electric potential, 



0(0) ~ v 



/in dz 



B T 

fj, l 



where I denotes the width of the electric current structures and L the length over which the integration 
should be carried out. Combining these estimates, and comparing the plasma velocity to the Alfvcn 
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velocity, we find 




) 



r) L 1 _ XL 
Hol d va S d 



where S is the Lundquist number based on the length-scale, I, of the current structures. 

Taking E n as an example, and taking L fixed (since in a realistic situation the length of a coronal 
loop should remain roughly constant with braiding), an approximate relation for d is d ~ 0.77 x 1CP"/ 2 
(as discussed in Section H) and, accordingly (taking L ~ 48, as for E 3 ) 



Therefore for any finite value of S we can find a value n such that s»dj and the n-braid cannot be in 
force-free equilibrium. In other words, the ideal approximation (S — > oo) in which a force-free equilibrium 
may be found for arbitrary n must eventually break down, even if the current structures themselves are 
of a large-scale. 

Within the more general theory of the magnetic braiding of the solar corona (where L may be, say, 
roughly the length of a coronal loop and remain approximately fixed with increasing braiding), from 
these models we expect the width of integrated parallel current layers, d, to decrease exponentially as the 
coronal field becomes more and more braided by photospheric motions. As we have shown, this increase 
in complexity need not cause a decrease in the scale of the current structures themselves, i.e. thin current 
sheets may not necessarily form (van Ballegooijen, 1985). However, after a finite amount of braiding d 
will become so small that, no matter how small the resistivity, the equilibrium cannot be maintained. 
One possibility is that the loss of equilibrium will be via magnetic reconnection. In three dimensions the 
rate of reconnection is given by J E\\dl where the integral is taken along the field line where the quantity 
has its maximum value. In order to determine where and how magnetic reconnection is triggered we 
would need to take a full 3D MHD simulation of the domain. However, the small-scale structures found 
in |J"_7||dZ| will necessitate a very high spatial resolution for any such accurate MHD simulation. For 
sufficiently high n and low rj this is beyond the reach of currently available computing power. Whether 
for moderate n and comparatively high r] effects due to the existence of small length-scales in J*j|idZ| are 
already evident is the subject of a current investigation. 
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6 Discussion and Conclusions 



In this paper we have readdressed the magnetic braiding process (Parker, 1972) in view of the considera- 
tion that in three dimensions it is the integrated parallel electric fields along field lines that are important 
for reconnection. In resistive MHD these are related to the integrated parallel currents, not just the 
current structure itself. The key outcome of the investigation is a theory that as the coronal field is 
braided the integrated parallel current along field lines will develop smaller and smaller structures. As a 
result we conclude that relaxation to a smooth equilibrium is not possible for a magnetic field with an 
arbitrary degree of braiding; magnetic braiding by photospheric motions will, at some point, lead to a 
lack of equilibrium in the coronal loop. 

In Section |2.1| a particular magnetic field configuration is constructed, modelled on the Borromean 
rings. We refer to this field as E 3 . The initial configuration for E 3 is taken as a magnetic field between 
two parallel plates containing six isolated twisted regions in an otherwise vertical field. Each twisted 
region has only a small amount of twist and the overall twist of the configuration is zero. Thus the braid 
can be taken to model a (straightened-out) coronal loop. Its aspect ratio is such that, together with the 
amount of twist, it would likely appear straight in coronal conditions. 

A Lagrangian numerical scheme was taken to relax the field E 3 towards a force-free equilibrium with 
the field on the boundaries held fixed. In the relaxed state the maximum Lorentz force inside the volume 
is two orders of magnitude lower, at ~ 10~ 2 , than that initially present. With the technique presently 
employed, numerical difficulties prevent further asymptotic relaxation towards j x B = 0. This is believed 
to be a result of the implementation of second-order differencing on a highly deformed grid and is the 
subject of a current investigation. Nevertheless, in the final state the twist in the configuration is evenly 
distributed along field lines and the current, initially in the form of six isolated closed current regions, is 
in two smooth, twisted, layers extending throughout the volume. Accordingly, our approximation to a 
force-free state is a reasonable one that allows us to deduce the main results. 

The equilibrium state is described in Section [3] There is no evidence for small- scales developing in the 
current structure, lending support to those arguing against Parker (1972) - for example van Ballegooijen 
(1985) and Craig & Sneyd (2005) - who claim smooth force-free equilibria are possible under arbitrary 
boundary perturbations. However, in three dimensions it is the presence of an integrated electric field 
component parallel to the magnetic field that is crucial for magnetic reconnection. In resistive MHD this 
is related to the integrated parallel current structure along field lines (J j\\dl), rather than the current 
itself. This consideration motivated us to consider the nature of J jndl in the braided field E 3 , both for 
the initial and relaxed states. 
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The variation in J*jndZ| between field lines originating on the lower boundary of the domain is exam- 
ined in Section |3j The \J j\\dl\ structure is found to be nearly the same in both cases, i.e. approximately 
conserved in the relaxation process, and there is seen to be significant variation in the value of this 
quantity between field lines that intersect the lower boundary of the domain only a small distance apart. 
We refer to this characteristic as 'thin integrated parallel current layers'. During the relaxation process 
the width of these layers is preserved, whilst the peak values in | J j\\dl\ are slightly enhanced, so creating 
larger gradients in \J j\\dl\ between the layers. 

In Section [4] the variation in the width, d, of these integrated parallel current layers with complexity, 
or degree of braiding, of the magnetic field is considered through an examination of the fields E n for 
various n. Here each increase in n can be thought of as being the result of additional braiding of the 
coronal field via photospheric motions. With increasing n, d is found to decrease exponentially. More 
generally then, the width of the integrated parallel current layers within a configuration is expected to 
decrease with increasing braiding. 

In Section [5] arguments are given relating to the consequences of narrow \J j\\dl\ layers on the pos- 
sibility of finding a force- free equilibrium. It is concluded that once L/d ~ S (where L represents the 
length of the braided field lines and S the Lundquist number) - a state that is inevitable as the magnetic 
braiding by photospheric motions continues - a smooth equilibrium state can no longer be attained. The 
precise manifestation of this lack of equilibrium is yet to be determined but we consider it is likely to be 
via magnetic reconnection events. 



18 



Bibliography 

Craig, I. J. D. and Sneyd, A. D., ApJ, 311, 451 (1986) 

Craig, I. J. D. and Sneyd, A. D., ApJ, 357, 653 (1990) 

Craig, I. J. D. and Sneyd, A. D., Sol. Phys., 232, 41(2005) 

Galsgaard, K. and Nordlund, A&A, JGR, 101, 13445 (1996) 

Hesse, M. and Schindlcr, K., JGR, 93, 5559 (1988) 

Hornig, G. and Priest, E., Phys. Plasmas, 10, 2712 (2003) 

Longcopc, D. W. and Strauss, H. R., ApJ, 437, 851 (1994) 

Longcopc, D. W. and Sudan, R. N., ApJ, 437, 491 (1994) 

Mikic, Z. and Schnack, D. D. and van Hovcn, C, ApJ, 338, 1148 (1989) 

Parker, E. N., JGR, 62, 509 (1957) 

Parker, E. N., ApJ, 174, 499 (1972) 

Petschck, H. E., in 'The Physics of Solar Flares', ed., Hess, W. N., p. 425 (1964) 

Pontin, D. I. and Galsgaard, K. and Hornig, G. and Priest, E. R., Phys. Plasma, 12, 052307 (2005) 

Priest, E. R. and Heyvaerts, J. F. and Title, A. M., ApJ, 576, 533 (2002) 

Schindlcr, K. and Hesse, M. and Birn, J., JGR, 93, 5547 (1988) 

Solanki, S. K. and Inhcstcr, B. and Schiissler, M., Rep. Prog. Phys., 69, 563 (2006) 

Sweet, P. A., in 'Electromagnetic Phenomena in Cosmical Physics, ed. Lehnert, B., p. 123, (1958) 

Van Ballegooijen, A. A., ApJ, 298, 421 (1985) 

Zweibel, E. G. and Li, H.-S., ApJ, 312, 423 (1987) 



19 



